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1.0.  INTRODUCTION 


It  was  observed  that  a  bump  existed  in  the  pressure-time  curve  during  the 
burning  of  solid  rocket  propellant.  The  pressure  bump  is  likely  due  to  a  con¬ 
centration  of  ammonium  perchlorate  in  the  propellant,  and  thus  related  to  the 
filling  process  [1].  A  common  feature  of  all  propellant  compositions  giving  a 
bump  is  a  high  solids  fraction,  not  less  than  70%  by  weight  (mostly  ammonium 
perchlorate)  dispersed  in  a  binder.  Unfortunately,  the  current  understanding 
of  a  densely  packed  particle/fluid  system  is  rather  poor  and  no  simple 
explanation  is  satisfactory.  Furthermore,  non-Newtonian  behavior  is  expected 
during  the  filling  process  which  makes  the  theoretical  analysis  of  the 
practical  system  very  difficult. 

In  a  previous  report  [2],  preliminary  investigation  was  performed  by 
using  the  multiple  velocity  field  model.  This  is  a  continuum  approach  as 
opposed  to  a  discrete  particle  track.  It  was  demonstrated  in  the  preliminary 
investigation  that,  by  using  different  velocity  boundary  conditions  for  the 
two  fluids,  partial  component  separation  was  predicted.  No  slip  was  assumed 
for  the  liquid  binder  while  free  slip  was  assumed  for  the  ammonium  perchlorate 
particles,  since  particle  velocity  may  have  a  finite  value  at  the  wall.  Even 
though  the  results  are  encouraging,  it  was  pointed  out  that  while  the  multi¬ 
fluid  model  is  more  complete  than  the  diffusion  model,  there  are  more  unknown 
coefficients  associated  with  the  model,  and  these  coefficients  need  to  be 
determined  either  experimentally  and/or  theoretically. 

In  this  report,  the  multifluid  model  will  be  adopted.  Emphasis  will  be 
placed  on  identifying  important  parameters  (dimensional  or  dimensionless)  for 
the  two-fluid  model.  An  attempt  will  be  made  to  narrow  down  the  parameters 
range  and  to  interpret  the  results  of  the  multifluid  model  that  are  relevant 
to  the  solid  propellant  casting  process. 


2.0.  IMPORTANT  PARAMETERS 
The  momentum  equations  for  a  two-fluid  system  are 

It  v.(e^p^UiUi) 

=  -  6^  VP  +  v.T^^  +  e^p^fj  -  k^2  (Uj  -  U2)  (1) 

(02P2U2]  (92P2lJ2^2^ 

-  -  02  VP  +  ^*^2m  ^  ^2^2^2  ""  ^21  ^^2  ~ 

where  0^,  pp  Uj  refer  to  the  volume  fraction,  the  material  density,  and  the 
velocity  of  the  binder  respectively,  subscript  2  refers  to  the  particles 
(ammonium  perchlorate),  is  the  shear  stress  of  component  1  in  the  mixture, 
fi  is  the  field  force  per  unit  mass  on  component  1  and  may  include  that  due  to 
shear  lifting  effect,  P  is  pressure,  Kj2  is  the  interfacial  momentum  transfer 
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coefficient  between  the  two  phases,  etc.  The  parameters  and  properties  needed 
are  f,,  f2»  02*  ^^2,  and  K2^.  Additonal  key  parameters 
are  trie  particle  size  and  its  distriliution.  The  volume  fractions  are  related. 


01  +  02  =  1  (3) 

The  field  forces  (fi  and  f2)  will  be  treated  separately  and  are  not 
included  in  the  present  investigation.  The  shear  stress  terms  and  T2m) 

are  not  readily  determined  at  this  time  for  dense  suspensions.  In  the  numer¬ 
ical  calculations,  it  is  simply  assumed  that  each  of  the  two  fluids  behave  in 
a  Newtonian  fashion  and  the  viscosity  of  each  fluid  is  not  affected  by  the 
presence  of  the  other  fluid.  Even  with  this  simplified  assumption,  it  is 
still  difficult  to  determine  the  viscosity  of  the  second  fluid,  i.e.,  ammonium 
perchlorate  particles. 

The  interfacial  momentum  transfer  coefficients  (K12  and  K2ij  can  be 
expressed  as 


12  "  ®lPl''l2  ’ 

(4) 

21  =  02P2''21  ’ 

(5) 

where  ?i2  inverse  relaxation  time  for  momentum  transfer  from  phase  2  to 

phase  1,  etc.,  and 

Ki2  =  K21  (6) 


from  action  and  reaction.  For  the  motion  of  a  particle  suspended  in  a  fluid 
in  the  Stokes  regime  [3], 


F  - 

•^21  p'2  ’ 


(7) 


where  a  is  the  radius  of  the  particle,  p,  is  the  viscosity  of  the  fluid,  and 
P2  is  the  material  density  of  the  particle.  Combining  Eqs.  5,  6,  and  7,  the 
following  equation  is  obtained: 
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a2 


‘'I 


(8) 


If  a  dimensionless  interfacial  momentum  transfer  coefficient  is  defined  as 


D  = 


(9) 


then 
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(10) 


D  = 


9  a. 
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which  is  valid  only  for  a  single  particle  suspended  in  a  fluid  medium  in  the 
Stoke's  regime. 

For  a  cloud  of  particles  suspended  in  a  fluid  stream,  the  inverse 
relaxation  time  is  [3], 


75  [1  -  0^)  1.75  p^u^ 

■"zi  "  “772  7”  ^  ~2  7 


(11) 


where  u^  is  the  superficial  fluid  velocity  based  on  unobstructed  flow  area. 
In  the  interest  of  the  solid  propellant  system,  the  largest  particle  size  is 
in  the  order  of  a  =  100  pm.  The  viscosity  of  the  binder  (p^)  is  relatively 
large  (-200  Pa-s  versus  lO"-^  Pa-s  for  room  temperature  water).  Thus,  the 
first  term  on  the  right-hand  side  is  dominating  for  the  solid  propellant 
system.  Therefore,  Eq.  11  can  be  written 


F 


21 


75  02  p^ 

^  * 

2  8^  P2  a 


Combining  Eqs.  4,  5,  and  12, 


0  = 


(12) 


(13) 


Equations  10  and  13  indicate  that  the  dimensionless  interfacial  momentum 
transfer  coefficient,  which  is  a  function  of  volume  fraction  only,  is  the 
relevant  scaling  parameter  for  interfacial  momentum  transfer.  For  given 
volume  fraction  (0i  or  02),  Eqs.  10  and  13  can  be  used  to  estimate  the  order 
of  magnitude  of  0  for  single  particle  and  a  cloud  of  particles  suspended  in  a 
fluid  stream,  respectively. 

In  addition  to  the  dimensionless  interfacial  momentum  transfer  coeffi¬ 
cient  D,  the  viscosity  ratio,  defined  as 


is  another  important  parameter  since  the  velocity  and  volume  fraction  distri¬ 
butions  will  depend  strongly  on  the  relative  magnitude  of  the  two  viscosities 
in  the  two-fluid  model. 

The  particle  size  and  its  distribution  will  significantly  affect  the 
outcome  of  the  solid  propellant  casting  process.  The  two-fluid  model  is  not 
capable  of  dealing  with  a  system  with  a  number  of  particle  sizes.  However,  by 
using  certain  physical  models,  such  as  the  Probstein-Sengun  bimodal  model  [4], 
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it  is  possible  to  include  the  contribution  frnn  various  particle  sizes  in  a 
two-fluid  model,  and  therefore  avoid  going  to  a  multifluid  model  which  will 
only  introduce  more  unknowns,  with  current  status  of  knowledge  a  fluid  system 
with  multiple  particle  sizes. 


3.0.  EFFECTS  OF  INTERFACIAL  DRAG  AND  RELATIVE  VISCOSITY 

investigate  the  effects  of  interfacial  momentum  transfer  and  relative 
viscosity  of  two  fluids,  steady-state  calculations  were  performed  for  an 
annular  flow,  which  is  the  two-dimensional  axisymmetric  flow  between  two 
concentric  cylindrical  surfaces  employed  in  previous  investigations  [2].  The 
characteristics  of  the  geometry,  values  of  properties,  and  boundary  conditions 
are  given  in  Table  1. 


Table  1.  Properties  and  Characteristics  of  Annular  Flow 


Inner  Radius,  R^ 

1.905  cm 

Outer  Radius,  Rq 

2.540  cm 

Gap 

0.635  cm 

Inlet  Velocity  Uj  =  U2  =  Uq 

1 .0  cm/s 

Inlet  Volume  Fraction 

02  (binder) 

0.2242 

02  (particle) 

0.7758 

Material  density 

P2  (binder) 

920  kg/m^ 

P2  (particle) 

1950  kg/m' 

Velocity  Boundary  Condition 

Bi nder 

no  slip 

Particle 

free  slip 

An  average  particle  size  with  a  radius  of  100  urn  is  assumed.  Since  the 
average  diameter  (or  radius)  is  used,  the  particle  volume  fraction  can  exceed 
the  maximum  packing  fraction  for  particles  with  a  single  diameter.  The  binder 
viscosity  is  assumed  to  be  equal  to  200  Pa-s. 

Figures  1  through  8  show  the  fully-developed  velocity  and  volume  fraction 
distributions  at  the  exit  of  an  annular  duct  for  different  values  of  the 
dimensionless  interfacial  momentum  transfer  coefficient  0,  and  viscosity  ratio 
Z.  For  a  viscosity  ratio  of  Z  =  1,  Figs.  1  through  4  indicate  that  the 
effects  of  interfacial  parameter  0  begin  to  appear  when  it  is  decreased  from 
0.50  to  0.05  and  the  effects  become  appreciable  when  0  =  0.005.  Significant 
separation  occurred  when  D  =  0.005,  as  shown  in  Figs.  3  and  4.  When  the 
interfacia]  momentum  transfer  coefficient  is  large  (D  >  0.5),  there  is  very 
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little  separation  and  the  volume  fraction  of  either  the  binder  or  the  particle 
remains  nearly  constant  across  the  duct.  This  is  understandable  since  the  two 
fluids  will  move  as  a  single  fluid  if  the  interfacial  drag  D  becomes  large. 

Figures  5  through  8  show  the  effects  of  viscosity  ratio  (Z  =  \i2_l\i\]  on 
velocity  and  volume  fraction  distributions  for  D  =  0.005.  It  can  be  observed 
that  the  viscosity  ratio  affects  the  volume  fraction  distribution  somewhat, 
but  the  most  pronounced  effect  appears  in  the  particle  velocity  distribution 
as  shown  in  Fig.  6.  If  the  particle  viscosity  is  much  larger  than  the  binder 
viscosity  (for  example,  Z  =  8.70),  then  the  particles  move  in  a  manner  as  if 
the  binder  is  absent.  On  the  other  hand,  if  the  binder  viscosity  is  much 
larger  than  that  of  the  particles  (for  example,  Z  =  0.04),  the  particle  veloc¬ 
ity  distribution  becomes  closer  to  the  binder  velocity  distribution.  The 
binder  velocity  distribution  appears  to  be  less  sensitive  to  the  viscosity 
rati  0. 

Two  additional  interesting  phenomena  were  observed  from  the  numerical 
results.  If  the  interfacial  momentum  transfer  coefficient  is  large  (D  >  1)  , 
the  volume  fraction  distributions  for  both  the  binder  and  the  particles  are 
uniform  and  are  independent  of  the  value  of  the  viscosity  ratio  Z.  This  means 
that  there  is  a  limiting  value  of  0  (si),  above  which  there  is  no  separation 
and  below  which  separation  does  occur.  This  result  indicates  that  it  is 
extremely  important  to  determine  the  value  of  the  dimensionless  interfacial 
momentum  transfer  coefficient  for  the  solid  propellant  system.  The  authors 
are  not  aware  of  any  existing  theoretical/experimental  model  that  can  accu¬ 
rately  predict  the  value  of  D  for  a  densely  packed  system  such  as  the  solid 
propellant  system.  Further  experimental  and  theoretical  investigation  in  this 
direction  will  be  very  helpful. 

The  second  interesting  phenomenon  occurs  near  the  transition  point 
between  separation  and  no  separation  (D  =  1).  '^’igures  9  and  10  show  the 
binder  and  the  particle  volume  fraction  distributions  for  D  =  1  and  Z  =  10. 
At  such  a  high  value  of  0,  separation  is  minute,  but  finite.  This  can  be 
observed  from  the  M-shaped  or  W-shaped  profile  and  the  small  increment  in  the 
vertical  scale  shown  in  Fig.  9  or  10.  The  reason  for  the  existence  of  an  M- 
shaped  (or  W-shaped)  profile  is  not  clearly  understood  at  this  time. 


4.0.  DISCUSSIONS 

In  the  preliminary  investigation,  both  the  multi  fluid  model  and  the 
single  velocity  field  model  were  proposed  [2].  In  the  multifluid  model,  the 
transport  properties  Fj^o*  Hm*  etc.,  are  not  readily  computed  or  measured.  The 
closeness  of  the  phase"  velocities  suggest  that  the  mixture  velocity  may  be 
sufficiently  representative.  This  leads  to  the  development  of  a  single  veloc¬ 
ity  field  model  described  in  Ref.  2.  It  is  recognized  that  the  inverse 
relaxation  time  [^12]  ^nd  the  mixture  viscosity  still  have  to  be  determined 
experimentally.  However,  it  is  difficult  to  visualize  how  a  single  fluid 
model  can  produce  particle  separation  observed  in  the  solid  propellant  system. 

In  the  investigation  of  dense  slurry  rheology,  Probstein  and  Sengun  [4] 
recently  proposed  a  bimodal  model  to  describe  the  polymodal  behavior  of  coal- 
water  slurries.  Their  model  appears  to  have  some  merit  and  seems  relevant  to 
the  solid  propellant  casting  process.  In  the  bimodal  model,  the  following 
definitions  were  employed 
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e  -  ^  ^ 

V  V +  +  V 

c  f  £ 

e 

®ff  V, 

®c  ■  V,  +  V  +17" 

f  £  C 

where  0^  is  the  volume  fraction  of  the  suspended  solids,  6^^-  is  the  fine- 
filler  volume  fraction,  is  the  coarse  particle  volume  fraction,  is  the 
volume  of  the  fine  filler,  is  the  volume  of  the  coarse  particles,  and  is 
the  volume  of  the  fluid.  It  is  further  defined  that 


where 

n^P  =  net  relative  viscosity  =  Uc/wq* 

Hj-p  =  coarse  relative  viscosity  =  W(-/pf» 

n^P  =  fine  relative  viscosity  =  Uf/uQ* 

=  apparent  viscosity  of  the  suspension, 

Uf  =  apparent  viscosity  the  suspending  fluid  with  fine  particles 
present,  and 

Ug  =  viscosity  of  the  pure  liquid. 

The  basic  assumptions  in  the  Probstein-Sengun  model  are 

’^cr  “  ^(®c) 


and 


Hfr  =  ^'(Qff)  • 


They  further  assumed  that  it  is  the  small  volume  fraction  of  the  colloidal 
size  particles  that  imparts  to  the  suspension  most  of  its  rheological  char¬ 
acteristics  (shear  rate  dependent  viscosity).  The  large  particles  are 
essentially  unaware  of  the  presence  of  fine  particles,  but  rather  see  a 
stiffened  fluid  that  has  the  same  viscosity  and  density  as  the  suspension, 
with  the  contribution  of  the  large  particles  to  the  viscosity  increase  coming 
about  from  hydrodynamic  dissipation  (independent  of  shear  rate).  These 
assumptions  were  supported  by  experimental  data  at  low  shear  rate.  The 
demarkation  between  fine  and  coarse  volume  fractions  for  six  different  coal/ 
water  slurries  ranged  from  2  to  8  microns,  with  an  average  of  about  4.5 
microns,  even  though  the  particle  diameters  in  the  slurries  ranged  from  0.5  to 
300  microns.  Thus,  the  viscosity  of  a  polymodal  system  can  be  represented 
approximately  by  the  viscosity  of  a  bimodal  system.  This  greatly  simplifies 
the  analysis  and  therefore  is  the  more  practical  approach  than  the  multifluid 
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model.  The  composition  of  solid  propellant  contains  particles  and  powders  of 
various  sizes  and  is  extremely  difficult  to  handle.  The  bimodal  model  of 
Probstein  and  Sengun  appears  to  offer  a  reasonable  alternative.  Of  course, 
there  are  limitations  for  the  application  of  the  bimodal  model  as  described  in 
Ref.  4.  In  addition,  the  bimodal  model  does  not  address  the  problems  related 
to  particle/particle  interaction,  shear  lift,  and  particle/  wall  interaction. 


5.0.  SUMMARY 

( 1)  A  two-fluid  model  is  employed  to  investigate  the  velocity  and  volume 
fraction  distributions  in  an  annular  duct.  A  no  slip  boundary  condition 
was  assumed  for  the  fluid  (the  continuous  phase)  and  free  slip  was 
assumed  for  the  particles  (the  dispersed  phase)  at  the  walls.  Two 
important  dimensionless  parameters  were  identifed  and  parametric  studies 
were  carried  out  by  varying  these  parameters.  The  most  important 

parameter  is  the  dimensionless  interfacial  momentum  transfer  coefficient 
D  (=  second  parameter  is  the  viscosity  ratio  of  the 

two-fluid  Z  [=  parameters  can  affect  the  velocity  and 

volume  fraction  distributions  significantly. 

(2)  Numerical  results  indicate  that  for  relatively  small  values  of  D,  the 

dimensionless  interfacial  momentum  transfer  coefficient  (D  <  0.05), 
significant  separation  between  the  binder  and  the  particles  are 

observed.  The  smaller  the  0,  the  larger  the  degree  of  separation. 

Viscosity  ratio  Z  also  affects  the  volume  fraction  distribution,  but  its 
influence  is  somewhat  less  than  that  of  the  interfacial  momentum  transfer 
coefficient  0. 

(3)  The  degree  of  separation  decreases  with  increasing  values  of  the 
interfacial  momentum  transfer  coefficient.  When  D  is  sufficiently  high 
(>1),  separation  no  longer  exists.  This  is  understandable  since  when 
interfacial  momentum  transfer  is  dominating,  the  two  fluids  move  like  a 
single  fluid  and  consequently,  there  is  very  little  separation. 
Currently,  there  is  no  theoretical  or  experimental  model  that  can 
accurately  predict  the  value  of  the  interfacial  momentum  transfer 
coefficient  for  the  solid  propellant  system.  Further  investigation, 
either  theoretical  or  experimenhtal ,  in  this  direction  will  be  extremely 
helpful . 

(4)  An  interesting  phenomenon  is  observed  near  the  transition  between  separa¬ 

tion  and  no  separation  (D  =  1).  The  degree  of  separation  is  rather  small 
at  D  =  1,  but  the  particle  volume  fraction  distribution  exhibited  an  M- 
shaped  profile  which  is  not  present  for  D  much  smaller  than  1.  The 

reason  for  the  existence  of  such  an  M-shaped  profile  is  not  clearly 
understood  at  this  time. 

(5)  The  bimodal  model  proposed  by  Probstein  and  Sengun  may  be  advantageous 
.  relevant  to  the  solid  propellant  system  since  it  allows  the  use  of  a 

■  lodal  viscosity  model  to  represent  a  polymodal  system.  This  is  a  more 

■  ..ctical  approach  than  using  the  multifluid  model  since  the  latter  will 
onl^-  introduce  more  unknowns. 
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(6)  It  is  suggested  that  a  logical  next  step  is  to  incorporate  the  non- 
Newtonian  behavior  of  the  binder  into  the  two-velocity  numerical  model. 
This  can  be  accomplished  by  modifying  the  present  computer  program  to 
include  the  effects  of  shear-rate  dependent  viscosity  in  the  following 
form 


u  =  c 


3r 


where  c  is  an  empirical  constant,  is  the  local  shear  rate,  and  the 
exponent  n  accounts  for  various  rheorSgical  behaviors  of  the  fluid. 
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Fig.  1.  Fully-Developed  Velocity  Distributions  at  the 
Exit  of  an  Annular  Duct  (Z  =  1) 
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Fig.  2.  Fully-Developed  Particle  Velocity  Distribution  at  the 
Exit  of  an  Annular  Duct  (Z  =  1) 


10 


er  Volume  Fraction 


1.000  1.050  1.100  1.150  1.200  1.250  1.300 


R/Ri 

Fig.  3.  Fully-Developed  Binder  Volume  Fraction  Distribution 
at  the  Exit  of  an  Annular  Duct  (Z  =  1) 
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Fig.  5.  Fully-Developed  Binder  Velocity  Distribution 
at  the  Exit  of  an  Annular  Duct  (D  =  0.005) 
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Fig.  6.  Fully-Developed  Particle  Velocity  Distribution 
at  the  Exit  of  an  Annular  Duct  (D  =  0.005) 
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Fig.  8.  Fully-Developed  Particle  Volume  Fr 
at  the  Exit  of  an  Annular  Duct  (D 
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Fig.  9.  Particle  Volume  Fraction  Distribution  at  the 
Exit  of  an  Annular  Duct  with  D  =  1  and  Z  =  10 
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Fig.  10.  Binder  Volume  Fraction  Distribution  at  the 

Exit  of  an  Annular  Duct  with  D  =  1  and  Z  =  10 
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